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Abstract. We compare evolutionary models for close-in exoplanets coupling irradiation and evaporation due re- 
spectively to the thermal and high energy flux of the parent star with observations of recently discovered new 
transiting planets. The models provide an overall good agreement with observations, although at the very limit 
of the quoted error bars of OGLE-TR-10, depending on its age. Using the same general theory, we show that 
the three recently detected hot-Neptune planets (GJ436, p Cancri, p Ara) may originate from more massive gas 
giants which have undergone significant evaporation. We thus suggest that hot-Neptunes and hot-Jupiters may 
share the same origin and evolution history. Our scenario provides testable predictions in terms of the mass-radius 
relationships of these hot-Neptunes. 

Key words, planetary systems — stars: individual (GJ436, p Cancri, p Ara) 
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1. Introduction 

Determining the evolution and formation mechanism of 
extrasolar planets in very small orbits, one to two hun- 
dred times smaller than Jupiter's orbit, remains a chal- 
lenging puzzle. This puzzle recently became even more 
complex when radial velocity surveys of unprecedented 
accuracy discovered three extrasolar planets similar in 
mass to Neptune, i.e. 14 to 21 M®, with orbital periods 
P=2.6-9.5 days and separations a=0.028-0.09 AU (Santos 
et al. 2004, McArthur et al. 2004, Butler et al. 2004). 
An emerging short-period, Neptune-mass, exoplanet pop- 
ulation presents us with intriguing questions about their 
compositions and origins, with potential implications for 
the understanding of our own solar system. Because of 
the presence of Jupiter-mass planets at orbital separations 
even closer than that of these hot-Neptunes, people ruled 
out the possibility that the latter were born as more mas- 
sive giant planets which have lost a significant fraction of 
their envelope. The current belief following this intriguing 
discovery was then that these hot-Neptunes are of differ- 
ent nature and have different formation and evolutionary 
histories than their more massive counterparts (Santos et 
al. 2004, McArthur et al. 2004, Mazeh et al. 2004). 

In a recent paper, we explored the effects of evapora- 
tion on the evolution of short-period jovian planets, us- 
ing a consistent treatment of the interior structure and 
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the irradiated atmospheric structure (Barman et al. 2001, 
Chabrier et al. 2004, Baraffe et al. 2004). The idea of evap- 
oration is supported by the recent discovery of an extended 
atmosphere around the transiting exoplanet HD 209458b 
(Vidal-Madjar et al. 2003, Vidal-Madjar et al. 2004), 
demonstrating the occurrence of strong atmospheric evap- 
oration and mass loss for short-period irradiated planets. 
These models could successfully reproduce the observed 
radius of the transit OGLE-TR-56 (Chabrier et al. 2004) 
but the case of HD 209458b remains difficult to explain 
(Baraffe et al. 2003), raising some contention over the 
general applicability of current evolution theories to hot- 
Jupiters. 

In this Letter, we first show that our models including 
irradiation effects from the thermal and the high energy 
fluxes of the parent star are consistent with recent observa- 
tions of four other new transits, confirming the peculiarity 
of HD 209458b. Motivated by this success and by the re- 
cent discovery of Neptunian exoplanets, we have extended 
our calculations to the evolution of irradiated and evap- 
orating gaseous planets below the mass of Saturn, down 
to masses comparable to those of Uranus and Neptune. 
We show that, contrary to present belief, all of the recent 
neptunian objects, along with their short-period Jupiter- 
mass cousins, could have originated from more massive 
irradiated gas giants which suffered from evaporation in- 
duced by the parent star's high-energy radiation. We also 
present characteristic mass-radius relationships that bear 
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the imprint of irradiation and atmospheric evaporation 
on a planet's evolution. These predictions will be tested 
in the very near future by surveys capable of detecting 
Earth-sized transiting planets. 

2. Analysis of newly detected transit planets 

Since the discovery of the two first and well confirmed 
transits HD 209458b (Charbonneau et al. 2000) and 
OGLE-TR-56b (Konacki et al. 2003), five new candidates 
were found very recently, providing more stringent con- 
straints on the theory of irradiated planets. The properties 
of these new transits and of their parent star are summa- 
rized in Table 1, which also provides the incident flux F- lnc 
irradiating the planet, as defined in Baraffe et al. (2003) 1 . 
Coupling consistently irradiated atmosphere models cal- 
culated with the appropriate incident fluxes (Barman et 
al. 2005) and interior structures, we have calculated the 
evolution of planets with parameters (a, M , -Fine) char- 
acteristic of each of the new observed transits. The evo- 
lutionary models and their input physics are described in 
Baraffe et al. (2003) and references therein. Figure^shows 
a comparison between evolutionary models and observa- 
tions for all transiting planets. OGLE-TR-lOb and HD 
209458b are displayed in the same panel, since they have 
similar properties, in terms of planet mass and incident 
flux, and the same model is appropriate to describe their 
evolution. When the age of the system is undetermined, 
we adopt an arbitrary age ranging from 1 to 5 Gyr. 

Besides irradiation, we also take into account the ef- 
fect of evaporation of the planet gaseous content due to 
the stellar incident high energy flux, based on the hydro- 
dynamic model of atmospheric evaporation developed re- 
cently by Lammer et al. (2003) and applied by Baraffe et 
al. (2004) to the case of the two first transits HD 209458b 
and OGLE-TR-56b. The description of the evaporation 
model and the details of its coupling with the planet evo- 
lution can be found in the two aforementioned references. 
Note that the evaporation rate is based on heating by the 
age- dependent stellar XUV and Lyman-a radiation, cali- 
brated to observed age-luminosity L xuv {t) I L\> \{t) rela- 
tionships (Lammer et al. 2003, Ribas et al. 2005), and 
takes root on previous studies devoted to escaping at- 
mospheres of the Solar system terrestrial planets (Opik 
1963, Watson et al. 1981). Such hydrodynamic evapora- 
tion yields mass loss rates in agreement with the observa- 
tionally determined lower limit for HD 209458b: M > 10 10 
g s _1 « 10~ 13 Mj yr _1 , at the present epoch (Vidal- 
Madjar et al. 2003). 

Models including irradiation and evaporation are 
shown by the solid lines in Fig. Q with initial masses of 
the planets ranging between 0.9 and 2.7 Mj. These ini- 
tial conditions remain uncertain, because of the uncertain- 
ties in the evaporation model, but the results hold, qual- 

1 Fi nc — i (^-) 2 Fi,, where iZ* and are respectively the 
radius and the total flux of the parent star, and a the orbital 
separation. 



itativcly, with different evaporation parameters. Figure ^ 
shows that consistent irradiated models, including irradi- 
ation effects due to both thermal and high energy fluxes 
of the parent star, are in excellent agreement with the 
observed radii at the proper ages for all of the planets ex- 
cept HD 209458b and possibly OGLE-TR-lOb, depending 
on its age. Figure fj] also shows that models with evapo- 
ration (solid lines) are almost undistinguishable from the 
case without evaporation (dash-dotted lines). This illus- 
trates the concept of critical mass defined in Baraffe et 
al. (2004), below which a gaseous planet should evapo- 
rate entirely after a certain age, for a given incident flux. 
The evolution of a planet above this critical mass, how- 
ever, remains unaffected by evaporation up to several bil- 
lion years. All the progenitors of the transiting planets 
fall in this latter domain except, interestingly enough, HD 
209458b and OGLE-TR-lOb, which are just below this 
limit. The specific case of HD 209458b was already dis- 
cussed in Baraffe et al. (2003) and Chabrier et al. (2004). 
These authors demonstrated, among possible explanations 
for the particularly large radius, that a mechanism, such 
e.g. as the one suggested by Guillot and Showman (2002), 
which would dissipate a fraction as small as 0.1%-0.5% of 
the stellar incident flux at the planet internal adiabat level 
could reproduce the observations of both HD 209458b and 
OGLE-TR-56b. We have presently tested this hypothesis 
on OGLE-TR-111, which is the least massive and the least 
irradiated planet (see Table 1). Adding 0.1% of the inci- 
dent stellar energy to the corresponding irradiated model 
yields a radius larger by 7% than the upper limit of the 
observed radius. If a dissipative mechanism, still to be 
identified, is a common feature of all close-in planets, our 
test indicates that either it has to be rather inefficient in 
OGLE-TR-111 or this planet has a significant rocky or icy 
core, yielding a ~ 5% smaller radius than a purely gaseous 
object (Saumon et al. 1996). 

3. The case of hot-Neptunes 

Current transits do not provide constraints on the evapo- 
ration rates, since, in their mass range, their mass-radius 
relationship remains unaffected by the process of evapora- 
tion. However, the good agreement of present models with 
most of the currently observed transits gives us confidence 
on a scenario based on a global process of irradiation and 
evaporation. We have thus extended this scenario to ini- 
tial masses below the critical mass in order to explore the 
regime of sub-jovian mass planets. 

Figure [21 displays the evolution of hydrogen/helium 
planets with initial masses ranging from 0.5 to 1 Mj or- 
biting a parent star under conditions characteristic of the 
three newly discovered Neptunian planets, namely 0.028 
AU from an M2.5 star (GJ436, Butler et al. 2004), 0.038 
AU from a G8 star (^Cancri, McArthur et al. 2004) and 
0.09 AU from a G5 star (/xAra, Santos et al. 2004). The 
observed spectra and magnitudes of the parent stars, com- 
plemented by well established age-activity relations for 
similar type stars, imply that these systems must be a 
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Table 1. Properties of transiting extrasolar planets and of their parent star. The following quantities are displayed: 
the orbital separation, the mass and radius of the planet, the mass, effective temperature, spectral type and radius of 
the star, the incident flux (see text) and the age when determined. 



Object a(au) M p (Mj„ p ) R p (R Jup ) M* (M ) T e //* (K) Sp. type R*(R ) log(F in c) age (Gyr) 



OGLE-TR-132 a 


0.0306 


1.19±0.13 


1.13 ±0.08 


1.35 


6411 


F 


1.43 


9.35 


0-1.4 


OGLE-TR^ 6 


0.023 


1.45±0.23 


1.23 ±0.16 


1.04 


6000 


G 


1.10 


9.26 


3 ±1 


HD 209458 c 


0.046 


0.69±0.02 


1 42+ - 10 


1.06 


6000 


G 


1.18 


8.72 


4-7 


OGLE-TR-10 d 


0.042 


0.57±0.12 


1.24±0.09 


1 


5800 


G 


1 


8.60 




OGLE-TR-113 e 


0.023 


1.35±0.22 


i ns+ - 07 


0.77 


4750 


K 


0.765 


8.53 




TrES-l / 


0.039 


0.76±0.05 




0.89 


5250 


K 


0.83 


8.32 


2.5 ± 1.5 


OGLE-TR-111 9 


0.047 


0.53±0.11 


1 00+° 13 

i - uu 0.06 


0.82 


5070 


G-K 


0.85 


8.12 





"Moutou et al. 2005 ^Torres et al. 2004 c Cody & Sasselov 2002 d Konacki et al. 2005 e Bouchy et al. 2004 
/ Sozzetti et al. 2004 9 Pont et al. 2004 



few billion year old. After ~ 2-5 Gyr, these planets have 
lost more than 90% of their initial mass, about 0.5 to 
0.95 Jupiter mass of gas, and ultimately become Neptune- 
mass planets. Note that this evaporation mechanism holds 
whether or not the planet has migrated inward from larger 
orbits since its formation. Migration implies the presence 
of a disk and upper limits for disk lifetimes are ~ 10 Myr 
(Armitage et al. 2003). Thus, regardless of migration, the 
planet will spend the majority of its life in its final short- 
period orbit, evolving under the influence of irradiation 
and evaporation. 

Given the large amount of mass lost by the planet, 
one may wonder about the consequences on the plane- 
tary orbit. A tentative answer can be inferred from an 
analogy with the orbital evolution of comets. If matter 
escapes isotropically, the planetary orbit remains unaf- 
fected. If mass loss occurs anisotropically, however, the 
orbit can be affected. Like for comets, mass loss will take 
place on the irradiated side of the planet and thus will 
occur initially in the direction of the parent star, exerting 
a force in the opposite direction and pushing the planet 
away from the star. This may affect the evaporation pro- 
cess and eventually quench it, if the orbital separation a 
increases sufficiently. A quantitative estimate of this effect 
will depend on the fraction of anisotropy of the escaping 
flow and its velocity. Such a study is beyond the scope 
of the present paper. However, the existence of close-in 
planets on very short orbits (a < 0.05 AU) suggests that 
this effect is small. This question definitly deserves a more 
thorough analysis. 

Note that in our model the upper atmosphere, where 
most of the incident XUV flux is absorbed, always re- 
mains within the Roche lobe radius of the planet. Above 
this limit, ~ 1.3 i?j for p 1 Cancri and ^iAra, Roche lobe 
overflow will occur, leading also to complete escape of 
the planet gas envelope (Gu et al. 2003, Lecavelier et al. 
2004). The runaway of the evaporation process (Baraffe 
et al. 2004) leads to the catastrophic expansion of the 
planets, as indicated by the dotted lines in the bottom 



panel, and eventually to the evaporation of the entire 
hydrogen/helium content. The mass-radius relationships 
predicted by our model are totally different from the one 
characteristic of a non-irradiated, non-evaporating, H/He 
Neptune-mass planet which, to first approximation, would 
have a radius R ~ 0.6 Rj after 1 Gyr (Zapolsky & Salpeter 
1969). 

4. Discussion and conclusion 

Our evolutionary models, accounting for both irradia- 
tion and evaporation processes, show that most of the 
close-in exoplanets discovered today, including the emerg- 
ing Neptune-mass population, may originate from signifi- 
cantly more massive Jupiter-like gas giants which have ex- 
perienced drastic mass loss due to the stellar high-energy 
radiation. This result is in direct opposition to the cur- 
rent belief that the hot-Neptunes are of different nature 
than the hot-Jupitcrs. This statement was used in a recent 
paper by Mazeh et al. (2004) who report an interesting 
anti-correlation between mass and orbital period of tran- 
siting planets. In their analysis, they ignored the recently 
discovered hot-Neptunes because of their supposedly dif- 
ferent nature. We show here that this supposition is by no 
means obvious and that Jupiter-mass and Neptune-mass 
close-in planets may share the same origin. Including the 
Neptune-mass planets in the sample of all short period 
planets weakens the case for the anti-correlation of Mazeh 
et al. (2004). 

Our scenario is certainly affected by the large uncer- 
tainties in the evaporation model of Lammer et al. (2003), 
which is based on the assumption that the planet un- 
dergoes maximal energy-limited evaporation. This idea of 
energy-limited evaporation, meaning that the escape rate 
is essentially determined by the amount of the high en- 
ergy flux absorbed in the upper planetary atmosphere, was 
originally suggested by Opik (1963), and further applied 
by Watson et al. (1981) and many others (see references 
in Lammer et al. 2003) to terrestrial planet atmospheres. 
As recently shown by Yelle (2004), energy-limited escape 
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log t 

Fig. 1. Evolution of the radius (in Rj) for currently known 
transiting planets. The incident stellar flux contribution, 
Fi nc , increases from top to bottom. The measured mass 
and orbital separation are indicated in each panel. The 
curves correspond to evolutionary models for three cases: 
without irradiation (dashed lines), with irradiation effects 
on the thermal structure but no evaporation (dash-dotted 
lines), and with irradiation and evaporation effects (solid 
lines). For the latter cases, the initial mass M, is also in- 
dicated in each panel. Note that in most cases the dash- 
dotted (no evaporation) and solid (evaporation) lines are 
undistinguishable (see text). 

rates for close-in giant planets may be limited by efficient 
cooling due to chemical species, such as H^" and H + , which 
decreases the efficiency of local heating due to the ab- 
sorption of the XUV flux by the upper atmosphere. The 
absence of detailed chemistry in the work by Lammer et 
al. (2003) may result in overestimated rates, as suggested 
by Yelle (2004). Indeed, this latter finds escape rates 20 
times smaller than the ones estimated by Lammer et al. 
(2003). With such significantly smaller rates, hot-neptunes 
can not originate from Jupiter- or even Saturn-mass plan- 
ets. Test calculations with an escape rate 10 times smaller 
than the one predicted by Lammer et al. (2003) indicate 
that initial masses less than ~ 0.1 Mj are required to 
reach a Neptune-mass with the orbital parameters of the 
three discovered hot-Neptune systems within a few Gyr . 




7 8 9 10 

log t 



Fig. 2. Evolution of evaporating giant planets reaching 
a Neptune-mass (indicated in the top panel by the long- 
dashed line, in Jupiter-mass units Mj) within a few Gyr. 
Three planet-star cases are shown for: the G5 star /i 
Ara (Santos et al. 2004) (solid), the G8 star p 1 Cancri 
(McArthur et al. 2004) (dashed), and the M2.5 star GJ 
436 (Butler et al. 2004) (dash-dotted). The planet's or- 
bital separation a (in AU) for each system is indicated in 
the upper panel. 

The formation of such low mass gaseous planets at such 
small orbital distances seems to be very unlikely within the 
framework of current formation models of giant planets. 
This certainly would cast doubt on our present sugges- 
tion of a common origin for close-in neptunian and jovian 
planets. 

The work by Lammer et al. (2003) and Yelle (2004) are 
the first attempts to understand all the complex processes 
involved in the evaporation of close-in giant planets. It is 
presently impossible to favor one of these models over the 
other. Since the Lammer et al. (2003) model yields maxi- 
mal energy-limited escape, it is important to examine its 
consequences on the planet fate as an upper limit case 
for the evaporation. The next step is to explore system- 
atically the effect of lower evaporation rates, as well as 
of other uncertainties inherent to our calculations, mainly 
the effect of initial conditions, of a rocky/icy core and of 
non-standard chemical composition in the envelope. This 
work is under progress. 

The discovery of neptunian-mass transiting planets 
will be crucial for testing the present scenario. As shown in 
§3, the proposed irradiation and evaporation global pro- 
cess will be directly revealed by the planet mass-radius 
relationship. As mentioned previously, a purely gaseous 
(H/He) Neptune-mass planet not affected by irradiation 
and evaporation during its evolution will have a radius 
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R ~ 0.6 i?j at present epoch (> 1 Gyr). Planets com- 
posed dominantly of solid material, like terrestrial plan- 
ets or the ice giants, Neptune and Uranus, will have sig- 
nificantly more compact structures (Neptune's radius is 
i~ 0.35 i?j at its present age). In contrast, if irradiation 
plus evaporation is a dominant process in the evolution 
of short-period exoplanets, with Neptune-mass planets or 
even smaller objects originating from Jupiter-like gas gi- 
ants, we predict significantly larger radii for a given mass, 
namely R > 0.8 Rj for the presently detected Neptune- 
mass objects, as illustrated in Figure |2 

The large radii we predict for Neptune mass planets 
place them well within the detection limit of most transit 
surveys and the discovery of such an object would strongly 
support our scenario. However, given the low probabil- 
ity of transit detections, no Neptune-mass transit discov- 
ery remains consistent with the low number of transits 
detected up to now. If no object with the presently pre- 
dicted mass-radius relationships is found, once larger tran- 
sit statistics become available, then several possibilities ex- 
ist. Perhaps the evaporation mechanism is not a dominant 
process in planet history, questioning the observed evap- 
oration of HD 209458b and/or the validity of the present 
evaporation rates. Or, conceivably, all the progenitors of 
Neptune-mass objects lie below the predicted critical mass 
(see bottom panel of Figure[2l and, at most, a rocky core 
remains after evaporation. This latter case would suggest 
that many of the short-period low-mass planets should be 
made mostly of rocks. 

On the other hand, transit detections of Neptune-mass 
objects with Jupiter-like radii would superbly corroborate 
the present scenario of irradiated and evaporating gaseous 
giants. Our ability to find transiting planets will improve 
significantly with the soon to be launched COROT and 
KEPLER missions. Comparisons between these observa- 
tions and the present calculations will potentially bring 
information about the composition, evolution, and the ori- 
gin of short-period exoplanets, unveiling part of their mys- 
tery and possibly changing the understanding of our own 
Solar system. 
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